Bidirectional Single-Electron Counting and the Fluctuation Theorem 
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We investigate theoretically and experimentally the full counting statistics of bidirectional single- 
electron tunneling through a double quantum dot in a GaAs/GaAlAs heterostructure and compare 
with predictions of the fluctuation theorem (FT) for Markovian stochastic processes. We observe 
that the quantum point contact electrometer used to study the transport induces nonequilibrium 
shot noise and dot-level fluctuations and strongly modifies the tunneling statistics. As a result, the 
FT appears to be violated. We show that it is satisfied if the back-action of the electrometer is 
taken into account, and we provide a quantitative estimate of this effect. 
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According to the second law of thermodynamics, the 
entropy of a macroscopic system driven out of equilib- 
rium increases with time until equilibrium is reached. 
Thus the dynamics of such a system is irreversible. In 
contrast, for a mesoscopic system performing a random 
trajectory in phase space and measured during a suffi- 
ciently short time, the entropy may either increase or 
decrease. The 'Fluctuation Theorem' (FT), which re- 
lies only on the microreversiblity of the underlying equa- 
tion of motion, states that the probability distribution 
P T (AS) for processes increasing or decreasing the en- 
tropy by ±AiS during a time interval r obeys the relation 



P T (AS)/P T (-AS) = exp(AS). 



(1) 



Remarkably, the FT remains valid even far from equilib- 
rium. It has been proven for thermostated Hamiltonian 
systems Markovian stochasticprocesses 0, 0], and 
mesoscopic conductors S S 0, Q ■ The FT is funda- 
mentally important for transport theory, one of its con- 
sequences being the Jarzynski equality [SJ, [l(| , which, in 
turn, leads to the 2nd law of thermodynamics. It also 
leads to the fluctuation-dissipation theorem and Onsager 
symmetry relations lll| , as well as to their extensions to 
nonlinear transport U 1, @, 0, II] • 

In electron transport experiments the entropy produc- 
tion is related to Joule heating, AS = qeVs/T, where 
q is the number of electrons (with charge e) transfered 
through the conductor during time t, and T is the tem- 
perature. Hence the FT can be formulated in terms of 
the distribution of transfered charge P T (q) at sufficiently 
long times, r > e/J, as follows 



Pr(q)/Pr(-q)^exp(qeVs/T) 



(2) 



The FT has been experimentally verified in chemical 
physics at room temperature [13] , while tests of the FT in 



mesoscopic transport experiments in milli-Kelvin regime 
have been lacking so far. On the other hand, recent ad- 
vances in time-resolved char ge d etection by a quantum 
point contact (QPC) [U, Hill [HI made it possible to 
measure the distribution P T (q) for single-electron tunnel- 
ing through quantum dots. This opens the possibility of 
testing the FT in mesoscopic transport. 

In this article we report on experiments on the 
direction-resolved full counting statistics (FCS) of single- 
electron tunneling in a double quantum dot (DQD) sys- 
tem, which is probed by an asymmetrically coupled 
QPC [3]. We analyze the experimental data in the frame 
of the FT and find that the form ([2|) appears to be vi- 
olated. However, it is still satisfied, if we replace the 
temperature by an enhanced value T* . We attribute the 
apparent overheating to the back action of the QPC de- 
tector and provide quantitative estimates of the effect. 
Moreover, we note that the form §\§ of the FT, valid 
for Markovian stochastic processes defined by general 
rates @, [j|, is satisfied, and we relate the entropy change 
to a ratio of the relevant nonequilibrium tunneling rates. 

Experimental test of the FT - The setup of our exper- 
iment [3] is shown in Fig. [TJ It consists of the DQD 
coupled to the QPC detector (Fig. QJi). The left and 
right gate voltages Vql , Vqr applied to the quantum dots 
are tuned in such a way that only three charge states of 
the DQD are allowed: \L), \R), and \D), denoting states 
where the left or right dot is occupied with a single elec- 
tron, or where both dots are occupied, respectively. Ac- 
cordingly, the current through the QPC, which is coupled 
asymmetrically to the DQD, switches between three dif- 
ferent values (Fig. [lb). This setup allows distinguishing 
electron tunneling in different directions. For example, 
the switching \L) — > \R) corresponds to the transfer of 
one electron from the left dot to the right one, while 
\R) — > \L) signals a transfer in opposite direction. The 
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FIG. 1: (a) Schematics of the system. Only one energy 
level in each dot, with the energies el and er respectively, 
contributes to the transport, (b) The QPC current switches 
between three values corresponding to the three charge states, 
(c) Test of FT ([2]) at several times. Lines with symbols: loga- 
rithm of lhs of Eq. © at several times; dashed line: q eVs /T* 
with T* — 1.37 K; dot-dashed line: qeVs/T. Inset: the dis- 
tribution P T (q) at t = 4 ms. 



time trace of the current taken over a sufficiently long 
time r allows one to determine the distribution of trans- 
fered charges P T (q). An example of such a distribution 
is shown in the inset of Fig. [TJ:. 

First, we perform a direct test of the FT, see Fig. [IF,. 
The combination \n\P T {q) / P T {— q)] indeed is a linear 
function of the transfered charge q. However, the slope 
approaches at long time r the value eVs/T* where Vs — 
300 /iV is the value of the applied DQD bias voltage, but 
the effective temperature T* = 1.37 K is a fit parameter 
(dashed line), which strongly exceeds the bath tempera- 
ture of the leads of T = 130 mK (dot-dashed line). 

To further test the time dependence contained in Eq. 
(J2) , we check the integrated form of the FT (Fig. \2k) , 



-qeVs/T 



(3) 
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Plotting both sides of Eq. ([3]) with the value of the elec- 
tron temperature, T — 130 mK, we observe a clear dis- 



2 4 6 
x (ms) 

FIG. 2: (a) Test of the integrated FT Q in the time in- 
terval < r < 8 ms. The time r is multiplied by the fre- 
quency of electron tunneling I/e ~ 370 Hz. Squares: lhs of 
Eq. © (denoted as P T {q < 0)/P T (q > 0)); dashed line: 
rhs of Eq. Q (denoted as (exp(— qeVs/T)) q >o); solid line: 
rhs of Eq. © with T replaced by T* = 1.37 K. All three 
curves are multiplied by exp(7Y/e) for clarity, (b) Normal- 
ized second, C2 = ((<? — (q)) 2 )/(q), and third cumulants, 
— ((<? — {l) ) 3 ) /{q)i °f the charge distribution P T (q). (c) 
The six transitions with F^ between three charge states. 



crepancy between both (open squires and dashed line). 
However, by adjusting the temperature in the right hand 
side of Eq. ([3]) to the effective temperature (solid line), 
we get a good fit for r > e/J. As we will discuss below, 
the apparent heating is caused by the back-action of the 
QPC electrometer, which operates in highly nonequilib- 
rium conditions. 

The experimental distribution of the transfered charge 
(inset of Fig. [TJ;) deviates strongly from a Gaussian 
shape. For illustration the second and the third cumu- 
lants of the charge distribution are plotted in Fig. [3t>- 
The normalized third cumulant remains close to 0.1 at 
all times. Since the FT §2§ is satisfied for any Gaussian 
distribution P T (q) , the existence of higher cumulants is a 
further indication of non-trivial behavior. 

Violation and recovery of the FT- Why does the FT 
appear to be violated? Quite generally, the FT for the 
system with four leads, depicted in Fig. [TJi, should be 
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formulated in terms of the joint probability distribution 
P T (q,q') of two charges eq and eq' transmitted through 
the DQD and the QPC, respectively 0], 

P T (q,q') = exp [(qeV s + q'eV QFC )/T] P T (-q,-q% (4) 

where T is the temperature of the leads. Since only 
the charge eq is measured, Eq. (j4|) should be integrated 
over eq' and afterwards, the right hand side reduces to 
exp[eVs/?1 Pr(—q) only if Vqpc = 0, which leads to the 
apparent violation of the FT. 

Then why is the FT recovered by introducing the ef- 
fective temperature T*? In the experiment, the DQD is 
in the sequential tunneling regime. Then the probability 
distribution P T (q) can be derived from a master equa- 
tion approach of the FCS 

SIS El- It be gins with the 

'modified' master equation 



d tP (t) = T(X)p(t), p T = {pl,Pr,Pd) 



(5) 



where p s is the occupation probabilities of the charge 
states (s = L,R,D). The transition matrix has a mod- 
ified form r(A) depending on the counting field A mea- 
suring the electron transfer through the barrier between 
the two quantum dots 



r(A)= 



T_RL — ^DL 
1 RL 6 



r M e- lA 

Tl_r — Tor 
Fdr 



Fld — F_R£) 



(6) 



The characteristic function (CF) takes the form [l7| 

ZtW =J2 a P T (q)e i « x = e T e rT Wp st , (7) 



where e T = (1,1,1) and p st is the stationary state, which 
is found from the equation T(0)p st = 0. In the long- 
time limit, T>e/J, the CF acquires an exponential form 
Z T {\) « e T:F ( x \ where ^(A) is the eigenvalue of the ma- 
trix T(A) with the largest real part. It satisfies 



= dct[r(A) - J - 1] = T z + KT 2 



drXrl^ld 



(e 1 ' 



D-r 



Dli LRX RD 



K'T 

r WD (c- iA -i) 



(8) 



where I is the identity matrix. Since the parameters K 
and K' are independent of A, we observe without solving 
Eq. ([SJ) , that the CF satisfies the identity, 

Z{\) = Z(-\ + ieV s /T*), (9) 
T* = eV s /\n[(T D RT RL T LD )/{T DL T L RTR D )].{lQ) 

Performing the inverse Fourier transform of Eq. ([9]) , we 
arrive at the FT in the form © with T replaced by T* . 
The argument of the logarithm in Eq. (fT0|) is the ratio 
of products of the tunneling rates corresponding to for- 
ward and backward cycles (i.e., the counterclockwise and 
clockwise cycles shown in Fig. [2]:). The logarithm of this 
ratio gives the entropy production associated with the 
transfer of one electron through the system 0, Q . While 



in the present example we find this transparent result we 
note that in more general systems with many transport 
cycles, it is in general not possible to define a unique 
effective temperature 

In the experiments the tunneling rates are estimated as 
T DR = 4kHz, Y RD = 0.3kHz, r DL = lkHz, T LD = 1.5kHz, 
Tlr = 1.7kHz, and Vrl — 1.8kHz. The values lie in 
a regime where they do not suffer from the finite band 
width of the QPC detector (~ 10 kHz). The effective 
temperature derived from Eq. (fTU)) . T* = 1.14 K, agrees 
well with the value 1.37 K obtained directly from the FT. 

QPC back-action - The tunnel rates are influenced by 
environmental effects. In particular, the DQD is neces- 
sarily coupled to the QPC, as modeled by capacitors Cl 
and Cr in Fig. [T^l. The nonequilibrium current noise of 
the QPC [l9[ produces fluctuations of the potentials of 
the quantum dots eSVr,,R, which in turn influence the 
tunnel rates 
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201 ] as known from the so-called P{E)- 
theory [2l| . Introducing three such functions Pj (E) with 
j = L, R, dd we find 

T RL = 27r\t dd \ 2 P dd {E L -E R ), (11) 
T DL = T R J dE f(E D -E L - vr - E)P R (E), (12) 

T DR = T L JdEf{E D -E R -Li L +E)P L (E). (13) 

Here t dd is the matrix clement describing tunneling be- 
tween the two dots, the rates Tl and Tr characterize 
the coupling between the dots and the leads. f(E) is 
the Fermi function, where we fix the chemical potentials 
of the leads as /.il = —Hr = eVs/2. The energies of the 
charge states El.rjj include the electrostatic energy of 
the electric field stored in the capacitors. The rates Tlr, 
r and Tro are given by the same expressions where 
the Fermi function / should be replaced by 1 — / and 
the argument of the functions Pj should be taken with 
opposite sign. 

The spectral functions PL/R/dd are expressed in terms 
of the phase operators <PL/R.(t) = /* dt' 8Vl / r{^') an d 
<fidd = <i>R-<pL as follows 



PAE) = 



x exp 



dt 



e iEt / e i<?j(t) e — i&/(0) 



dt 



AEt 



2?r \ / J 2tt 

d^e^So.jH + S^HKe-^-l) 
2tt lo 2 



(14) 



Here Soj(oj) is the spectral density of thermal fluc- 
tuations of equilibrium environments, including the 
impedance of the external circuit Z cxt (Fig. UK), 
etc. The part Sy^ within the Gaussian approxima- 
tion is proportional to the non-equilibrium and non- 
symmetrized current noise of the QPC as S$jJ C (w) = 



4 



K]\Z t {u)\^ pc {u) 



where 



1 7T 



T(l-T)(cj±eV Q p C ) , 2T 2 



, ± 



1 _ e -(w±eV Q P C )/T 



+ 



1 - e-^/ T 



Here T is the QPC transparency and Z t is the impedance 
of the electromagnetic environment seen by the QPC; 
Z t {u>) = l/(— iu)C + 1/R) where the resistance and ca- 
pacitance are R=l/ '(-^Qpc+^cxt) and C= {3C +C G )/2+ 
ClCr/{Cl + Cr). The factors Kj characterize the cou- 
pling between QPC and quantum dots given by certain 
ratios of the capacitances. 

When the QPC is in equilibrium Vqpc = 0, the func- 



tions Pj obey the detailed balance Pj(E)/Pj 



etc. 



E) = 
From 



e E / T and thus T RD /T DR 
Eqs. (|10M13p it follows that in this case T*=T. However, 
in the experiment, the QPC is biased at a rather high 
voltage e Vqpc 3> T and the tunnel rates through the cen- 
tral barrier are approximately T RL ^ LR ^T m&K /[l + (E R — 
E L f/TH where r = T(l - T)(e 2 K dd R) 2 eV QPC /2n 
and r max = 2 |id<j| 2 /ro- From the experimental values 
r ma x ~ 7 kHz, T « 30AieV, Vqpc = 800 fN and T= 0.19, 
the parameters can be roughly estimated as \tdd\ ~ 30 
MHz, KddR ~ 5 kO. The latter value is in agreement 
with other experiments, see i.e. [151 ] . Since the difference 
between the tunnel rates in opposite directions is signifi- 
cantly reduced, Tl R ^^ R l, we observe that the effective 
temperature T* is enhanced. From Eq. (|10[) we roughly 
estimate T* ~ eVs/ ln^e^/To] ~ 1 K, in good agree- 
ment with our previous findings, thus further supporting 
the QPC back-action model. In order to reduce T* one 
can reduce Vqpc or the external impedance seen by the 
QPC R. 

Other environmental effects do not change our find- 
ings by much. In GaAs nanostructures acoustic phonons 
modify the tunneling properties via piezoelectric and de- 
formation coupling. At experimental values of the QPC 
current /ape ~ 12nA the phonons stay almost in equi- 
librium 16[ . Therefore, the phonon effect is absorbed in 
the equilibrium part of the spectral density Soj, and does 
not affect the FT. An intrinsic back-action of the QPC 
is often addressed in the context of quantum measure- 
ment 22(. However, in the present case we estimate this 
effect to be negligible. Detailed discussions of these envi- 
ronmental effects and derivations of the theory presented 
in this article, based on the real-time diagrammatic tech- 
nique [jjj UK , will be published elsewhere [13] • 

Summary - We have investigated experimentally and 
theoretically the fluctuation theorem for Markovian 
stochastic processes by studying the direction-resolved 
full counting statistics in a double quantum dot system 
via a nearby quantum point contact electrometer. We 
found an apparent violation of the fluctuation theorem, 
which we attribute to the non-equilibrium electromag- 
netic fluctuations generated by the shot noise of the quan- 
tum point-contact electrometer. We also demonstrated 



that the FT is recovered if we adopt an effective value 
for the temperature ten times higher than the electron 
temperature. This effective temperature depends on the 
entropy production associated with the transfer of one 
electron through the double quantum dot, which in turn 
can be expressed by a ratio of forward to backward tun- 
neling rates. 
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